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Cm Abstract. We introduce a simple method for the investigation of the classical stability of static 

^ solutions with a horizon in Lovelock gravity. The method is applicable to the investigation of 

Q high angular momentum instabilities, similar to those found by Dotti and Gleiser for Gauss- 

'^^ Bonnet black holes. The method does not require the knowledge of the explicit analytic form of 

the black hole solution. In this paper we apply our method to a case where the explicit solution 

^^ is known and show that it identifies correctly the resulting unstable modes. 



Q^ 1. Introduction 






X 



^ 



Theories formulated in more than 4 space-time dimensions have been used to provide a unified 

description of particle interactions. Kaluza [1] and Klein [2], in their celebrated papers, unified 
^^ electromagnetism with gravity in four dimensions by starting from five dimensional Einstein 

^ action and compactifying one space dimension on a circle. Schwarzschild black hole solutions 

lO were generalized to D dimensional Einstein gravity by Tangherlini [Jij. Once it was realized 

C^ that superstring theories, candidate quantum theories of gravity capable of unifying all known 

L^ forces, are consistent only in ten dimensions the study of higher dimensional models became 

more important and relevant. Ten dimensional Einstein action arises in the zero slope limit of 

(<— ^ superstring theory. However, perturbative contributions to string theory generate corrections 

T-H to the Einstein action [5] . The correction terms contain higher powers of the curvature tensor. 

^~~^ Now most gravity theories with higher order curvature corrections have ghost excitations, with 

^ the notable counterexample, the Lovelock Lagrangian 0. The equations of motion in Lovelock 

gravity contain only second order derivatives and the metric is consistent with unitarity. We 

can write the Lovelock Lagrangian as 



n=0 

where the terms, Ln, are constructed from antisymmetric combinations of powers of the 

curvature tensor as 

1 " 

T _ ^j:aibi,a2b2,...,a„b„ TT p adi f<y\ 

4=1 

where Rab'^ is the curvature tensor and the generalized Kronecker delta symbol. 6, is the unique 
tensor that is antisymmetric under the exchange of pairs atbi and ajbj or pairs of Cidi and Cjdj 
and symmetric in the exchange of any ai and bi or a and di. The n = term of M corresponds 
to the cosmological constant, while the n = 1 term is Einstein gravity (qi = 1, and Li = R). 



The term L2 is the Gauss-Bonnet term. Only a finite number of terms, [{D — l)/2], contribute 
to the equation of motion at any fixed D. 

As a rule, actions are built based on general principles, such as symmetry and causality. All 
terms satisfying these requirements and built from the appropriate fields should be included in 
the Lagrangian. In this sense, there is no reason, why higher order Lovelock terms should be 
excluded from ([I]). As the dimension of the ratio of coupling constants an+i/an is [length^], 
higher order terms are of higher dimensions and, as a rule, their contributions become important 
at shorter distances. In other words, solutions in Lovelock gravity reduce to solutions in Einstein 
gravity when r^" >> a„ for all n < [{D — l)/2]. 

Analytic spherical black hole solutions in 5 dimensional Gauss-Bonnet (GB) gravity were 
found by Boulware and Deser [7J. The surprising property of GB black holes is the existence 
of a minimum mass Mmin, proportional to the GB coupling, 02- At the minimal mass the 
radius of horizon shrinks to zero and the Hawking temperature vanishes, making minimal mass 
black holes thermodynamically stable [8j. However, Gleiser and Dotti have shown [9j, based on 
the decomposition of excitations by Kodama, Ishibashi, and Seto [10] |llj that below a critical 
mass, Mc > Mmin GB black holes have a classical instability against high harmonics excitations. 
To find the location of instability in j9] a Schrodinger-type differential equation was derived 
and solved for linearized excitations. The equation is diagonal in, and as usual, the potential 
is dependent on the five dimensional generalized angular momentum, L. Thus, Mc is also 
dependent on L. Beroiz, Dotti and Gleiser [12], pointed out that the onset of overall instability, 
which is attained at the maximum of Mc{L) over all admissible L values, is reached in the 
L — )• 00 limit. Konoplya and Zhidenko [13j showed that there are the no other instabilities 
of the 5-dimensional GB black hole than the ones found by Dotti and Geiser. In the present 
paper we apply a simplified method to the investigation of the classical stability of 7-dimensional 
spherical symmetric black holes in n = 3 Lovelock gravity. The exact solution for such black 
holes was found in [14] . Takahashi and Soda generalized the formulation of Geiser and Dotti 
and found a master equation and general conditions for the stability of higher order Lovelock 
theories |15J|16] for positive Lovelock couplings. Camanho, Edelstein, and Paulos completed a 
more general investigation of the stability of higher order Lovelock black holes |I7J . 

All the above investigations rely on finding an analytic solution to the Lovelock equations. 
However, other static solutions with an event horizon (dubbed 'black objects' below), for which 
such exact solutions are not known, are also of interest. The principal aim of this paper is to 
demonstrate that a simple analytic method for the investigation of high harmonics instabilities 
can be applied and high harmonics instabilities can be found even in the absence of analytic 
solutions of the Lovelock equations. We will briefly discuss an example for a geometry in which 
such a problem arises at the end of the next section. Using these examples, we wish to check the 
correctness of our conjecture that 'inert' black objects do not exist even in the limiting sense. 
Stated more precisely, suppose that the radius of horizon of a black object of a certain topology 
has a range < r/j < 00. If its AdM mass is finite and its Hawking temperature vanishes in 
the limit Vh — ?• 0, than it approaches an inert state. Lovelock black holes for positive Lovelock 
couplings do have these properties. Our conjecture is that for all black objects, when they exist, 
the inert state is always shielded by a classical high harmonics instability at a finite radius of 
horizon. In other words, before the object reaches its inert state by Hawking radiation it will 
undergo a catastrophic disintegration due to the classical instability. 

The classical stability of black holes and other black objects is not merely of academic interest. 
Sufficiently low scale gravity models could imply the production of black holes at future particle 
accelerators. If space-time is more than four dimensional, then Lovelock terms change the 
spectrum of observable black holes. In particular, classical instabilities at nonzero mass would 
significantly effect the experimental signature of black hole production. [18] 



2. Large angular momentum instability 

We will study the metric 

ds^ = gAB dx^ dx^ = -fir) dt^ + -^dr"^ + r^dnD-2 (3) 

representing a D dimensional static black hole. Linearized perturbations, qab — ^ SAB + ^AB; 
can be classified into scalar, vector, and tensor types [lO] [15] [16] . In the present paper we will 
concentrate on scalar perturbations. For scalar perturbations Hab = S r]ABii",t), where S* is a 
spherical harmonic on S^~'^. In other words, S satisfies the Laplace equation corresponding to 
metric g 

AgS = -PS = -L{L + D-3)S. (4) 

As the unperturbed metric is static, the coefficients of the equation of motion are independent 
of t and a solution for Hab, having exponential dependence on time can be found. In other 
words, we can write 

hAB = e^'fAB{r,g)S. (5) 

Gauge transformations allow us to bring /ab to a simple form, whereas the non-vanishing 
components of Jab are fu{r), ftr{r), frr{r), and fsir), where we defined fij{r,g) = gij fs{r) 
[9]. If the system of linearized equations for the perturbations admits a solution with (x)^ > 
(runaway solution) then the metric gAB is unstable. 

Geiser and Dotti derived a Schrodinger-like differential equation for a single component of 
Hab for the D = 5 GB black hole metric. The energy eigenvalue of the equation was E = — w^, 
where uj was defined in ([5|. As usual, the Schrodinger potential depended on L^. They found test 
functions, such that the expectation value of the Hamilton operator was negative, establishing 
the existence of runaway solutions. Such solution existed, provided the ADM mass of the black 
hole was smaller than a critical value, Mc > Mmin- At the critical mass the instability was 
restricted to the horizon. If the mass was smaller the instability was extended up to a radius, 
outside of the horizon. The critical mass depended on angular momentum, taking its largest 
value at L — t- oo. This fact was used in a subsequent paper, [12], to find a simplified method of 
finding high angular momentum instabilities, which only relies on the limiting behavior of the 
potential at L — )• oo. After the elimination of all but one component of the perturbations the 
equation can be brought to the form 



^w^ 



f{r) — h{r) 



+ F/i(r) = -6j2/i(r), (6) 



where the potential has the following asymptotic behavior at large L |12j : 

V = L^f{r)vL{r) + vo{r) + ..., (7) 

Then it is easy to see that if we find an ro, such that /(ro)t'L(ro) < then the lowest eigenvalue 
of (pi) must be negative as w^ ~ O(L^) for large L. Then the static solution is unstable. Near 
the horizon /(r) ~ [r — rh)f'{rh) (where rh is the radius of horizon) and the question of stability 
hinges on the sign of f'{rh)vLi'i"h) ■ Then the onset of the instability is determined by the largest 
zero oi f'{rh)vL{rh)- 

A significant complication presents itself when applying the above method to solve the 
stability problem when the system of second order differential equations for perturbations cannot 
be decomposed to yield a single second order equation. For such a problem we propose using 
an alternative method of expanding the differential equations first in L~^ and f^~^, keeping 
the leading order contributions only, and then solving the truncated equations for a single 



component. This can provide an expression for wl(^), which ahows one to tackle the stabihty 
problem. We will apply this method to the investigation oi D = 7 black holes in Lovelock 
gravity though in this particular case we have been able to find the exact potential near the 
horizon. To prove our point we will avoid using the explicit form of the known solution during 
our investigation. We will only use properties derived from the equations of motion, such as 
asymptotic expansions and horizon expansions, from which one can determine the potential vl 
near the horizon. 

As an example, it is instructive to apply our method to the study of the classical instability 
of black branes in Lovelock gravity. These objects are defined as extensions of black holes into 
n additional compact dimensions with metric 

ds"^ = -f{r)dt^ + jfldr"^ + r^dQd + h{r)-fijdw'dw^ , (8) 

where jij is an n dimensional metric for compact coordinates, /(r) is assumed to have a zero at 
some r > 0, and h{r) is bounded. Though analytic solutions for black branes are not in general 
known in Lovelock gravity, a numerical investigation of the static solutions is quite feasible. The 
numerical investigation can provide a one-to-one correspondence between the radius of horizon 
and the ADM mass. Then one can perform an analytic investigation of the instability in the 
infinitesimal neighborhood of the horizon because the metric functions have a unique analytic 
horizon expansion. Thus, one can also determine the horizon expansion of VL{r) in arbitrary 
order of r — r^. This, combined with the numerically determined correspondence between the 
radius of horizon and the ADM mass allows one to find the onset of instability as a function of 
mass. We will investigate the instability of black branes in Lovelock gravity, using this method, 
in a subsequent publication. 

3. Instability of D = 7 Lovelock black holes 

At -D = 7 the first four terms of series M contribute to the equations of motion. An 5*5 
symmetric black hole solution depends on the two constants, 02 and 03, in addition to the 
Newton's constant, and the ADM mass. Static, spherically symmetric black hole solutions in 
7-dimensional Lovelock gravity were found by Dehghani and Shamirzaie [14j. The solutions of 
black holes with third order Lovelock terms were classified by Camanho, Edelstein and Paulos 

m- 

For D = 7 metric (JSJ) leads to the following equation of motion: 

S't = S: = A {fir) c{r) - b{r)]} = 0, (9) 



2^5 



where 



c{r) = r4 + 24r2a2[l-/(r)]+72a3[l-/(r)]^ 

b{r) = 4r[l-/(r)][r2 + 602.(1 -/(r)). (10) 

Depending on the choice of Lovelock couplings 02 and as there are at most three solutions 
of (§ 

2(6q3 + r^a2) - z (y) r\as - 2al) + z*{2 Vf'^ 



fir) 



1203 
where 



111 



Y = W+^W^ + 2r^^{a3-2al), (12) 

W = r^a2 i^al - 803) - 18a|m, (13) 



and the three possible solutions are defined by the choices z = 1, e^*'^/^, e 2 2 7r/3_ 
m is the coefficient of the asymptotic expansion series 

/(0^1 + «r-2-^ + A + ... (14) 

m is related to the ADM mass, MadMj as 

m = MADM^ — • (15) 

As m is proportional to Madm we will refer below to the former as "mass" , as well. The solution 



is asymptotically Minkowski if a = in (14), while it is AdS (dS) type if a is positive (negative). 



Then the radius of horizon (when there is a horizon exists) is 



r/i = \ —6a2 + v/36a2 — 24a3 + m, (16) 

Now we turn to the analysis of the stability of the D = 7 black hole in Lovelock gravity, 
using our simplified method. As we indicated above, we will refrain from using the explicit form 
of the black hole solution, given by (11)-(16) but rather rely only on ([9]) and expansions based 



on it. These are readily available even for metrics for which the equations cannot be solved. 
Let us first find the coefficients of the asymptotic series expansion in r. It is easy to see, by 
replacing /(r) in (|9|) by an asymptotic series in powers of r~^, that the only possible leading 



order nonzero coefficients are those given in (14). Furthermore, we find that the only possible 
solutions for coefficient a are 

'0, asymptotically Minkowski solution, 

^f^^-. L^2 _ 2 



4^ (q;2 — V'^2^ — a'^s) ' asymptotically AdS or dS solution, (17) 

4^ (02 + a/o2^ — f as) ; asymptotically AdS or dS solution. 

The constant m (the "mass") is undetermined, but once it is fixed it determines all further 
asymptotic expansion coefficients for a given choice of a (as long as that value is admissible for 
the given Lovelock couplings). Assuming a finite radius of convergence, and no naked singularity 
(for real r), a and m determine the solution uniquely. We will find that in a similar way the 
horizon expansion is also defined uniquely for given radius of horizon, r/i, and Lovelock couplings. 
\n\ shows that besides a solution with Minkowski asymptotics, which exists for all Lovelock 
couplings, one solutions with AdS and one with dS asymptotics exist if as < 0. Furthermore, 
if < as < 3 a2^ / 2 then two AdS solutions exist if 02 > and two dS solutions if 02 < 0. 
This result agrees with the classification of [20]. 

The horizon expansion of /(r) can also be found from ([9]), without the knowledge of an 
analytic solution. In particular, one obtains 

fn . _ K^ _ 4rfc(rg + 6a2) 



c{rh) 4 + 24r2a2 + 72as ' 

^^""^ - Hr,W ~ ^r^ + 24r2a2 + 72a3)3^'^'^ +'°^'^"^ 

+ 432r^(2a^ - as) - 15552r^a| - 31104a2a| + 864r^a2(8ai - 9as)] . (18) 

Higher order expansion terms can be readily obtained in a similar manner. The horizon 
expansion is unique for given Lovelock couplings and a given choice of r^. Again, if there 



is no naked singularity at or above the horizon the series has a finite radius of convergence and 
can be continued analyticahy to all positive r > r^. Consequently, it will coincide with the 
analytic continuation of one of the asymptotic solutions. Then, by necessity, there will be a 
one-to-one map of admissible values for r^ and of mass parameters, m > 0. This map cannot 
be determined analytically, without the knowledge of the exact solution, but can be obtained 
by numerical integration of the equation(s) of motion. 

As is seen from (17) if as < 802^ then there are three asymptotic solutions of (|9|) . However, 



for every choice of r^ there is an unique horizon expansion. Then two of the three asymptotic 
solutions do not have horizons. They must end in naked singularities at some finite r^ > 0. In 
case the equation of motion cannot be solved analytically numerical integration starting from 



the horizon selects which of the three asymptotic behaviors (17) corresponds to the solution 
with a horizon. 

A further comment, concerning the horizon expansion is that the expansion diverges at 

co{rh) = rt + 24rla2 + 72a3 = 0. (19) 

The largest root of this equation for r^"^ > is 

rs"^ = -12a2 + \/l44a2^-72a3. (20) 

We define co(r) = c(r) |j(r)=o- As we will see below, for a given choice of the Lovelock couplings 
rg provides a lower or upper bound on physical values of r^- We prove now the following (only 



to cases when the right hand side of (20) is real): If at all real r > r^ the solution is regular 
then c{r) does not change sign. The proof is very simple: If there is no singularity above the 
horizon then f'{r) is bounded at any finite r > r^, Then it follows from ([9]) that c{r) = implies 
5(r) = 0, as well. The only simultaneous solution of these two equation for r"^ and /(r) is r^ = 
and /(r) = 1. 

In particular c{rfi) > for a solution with Minkowski asymptotics, because then /(r) — t- 1 
and c{r) ~ r^ > when r — )• 00. Similarly, one can see that if as < any regular solutions with 
dS or AdS asymptotics must have c{r) < 0. In particular, c(r/j) < 0. However, as shown by 
( |19| ), for large r^ c{rh) > 0. Thus, at as < for asymptotically dS or AdS solutions the radius 
of horizon of a black hole is limited to r^ < rg. For Minkowski asymptotics the reverse is true, 
Ts > rh- 



Returning to (20), it follows that the r/^ — t- limit is possible only if as > 2a2^ (Minkowski 
asymptotics only) or < as < 2a2^ and a2 > (Minkowski or AdS asymptotics). Since our 
main interest is the fate of black holes (or black branes) in the limit of r/^ — t- we will concentrate 
on solutions at as > 0, though instabilities can also be analyzed in a similar way when as < . 

Writing the equation of motion, linearized in time dependent perturbations as 



^4^G^ + nl (21) 



the equations for the perturbations take the general form 

^A = E br,L,u. fAB + 7',L,o. JaB + 7.%,. Iab] = 0, (22) 



A,B 



where the coefficients % i^, k = 0,1, 2, depend on the metric components of the static solution. 
Rather than solving this system of equations exactly we truncate these equations, by expanding 
IrLu] ™ 1 / L'^ and in \ j uj^ and keeping only the leading order term of the expansion. The 



truncated equations have the form 
L2 



nl 



njr 
rL~. 



^^^^rf{r)c{r)frr + c'{r)fs] + ^-^^^^[f<^-rfir)ac^O, 



c{r) 



2r6/(r) 



[buj^fs-lQu:^rf{r)Ut-L'ftt] 



d{r) L 



2r 



T—fs 



{2 r c(r) /;, - [2c'(r)/(r) + c(r)/'(r)] /^} c^ 0, 



•1/6 



4 j-6 
/Wc'(r) 



/rr 



c'(r) „ ^ c"(r) ^ 



8r5/(r) 



(23) 



where c(r) was defined in (10). 

The rest of the equations, Ti^ = 0, are not independent from the above four. We present the 
equations without truncation in the Appendix. 



Now taking the combinations of equations (23) 



4cW F'2^ei 



bc[r) u)"^ j[r) rc[r) ^ 

we obtain the following second order differential equation for fs- 



fir. 



d 
dr 



miss 



+ £'nr)vL{r)fs = -•J'fs, 



where 



viir) 



2 [c'jr)]'^ - c{r)c"{r) 
5r c(r) ^{r) 



(24) 



(25) 



(26) 



Note that (26) coincides with the eq. (48) of [15j derived using the explicit solution of the 
Lovelock equations. A term, proportional to f'g, and finite in the limit L,u} ^- oo, has been 



omitted from (26) as it can be eliminated by factoring out a multiplier from fs- This would 



only contribute by a sub-leading, 0(1), term to the potential. 

Note that our method enabled us to derive the functional form of VL{r) without assuming 
the exact functional form of the Lovelock black hole solution. Furthermore, dropping the 0(1) 
derivative terms from (23), one obtains four algebraic equations for the perturbations /ab, of 
which Til = ^^^ ^r — imply that in leading order the perturbation frt vanishes. Then 
dropping frt, along with other non-leading order contributions from equations Tir = and 



K 



v> 



we obtain a system of two linear algebraic equations for fu and fs- Eliminating fu 



leads to the single equation for fs which reads as 



L^f{r)vL{r)fs = -oo^fs, 



(27) 



with VL{r) given in (26), with c{r) replaced by co(r). 



To summarize, our procedure consists of the following steps: (1) Derive the system of 
equations for linear perturbations form the Lovelock equations (see the appendix for the complete 
equations), (1) Drop all derivative terms from the system, (2) Expand the coefficients of 
perturbations fAB{f) to leading order in L~^ and w"^, (3) Discard equations in which the 
only leading order terms contain off diagonal perturbations, and omit the terms containing off 
diagonal perturbations from the rest of the equations, (4) Calculate the determinant, Det, for 



a subset of independent linear equations for the diagonal perturbations fAA^r), (5) Solve the 
equation Det = for w^, to obtain an equation for the potential 



u' = -f{r)L'vL{r). (28) 



Equation (27) verifies that our simplified procedure for finding vl provides the same leading 



order potential as [I3],[TB]). While this method of finding the potential does not constitute a 



significant simplification compared to the "exact" derivation of (25), it is the only way to derive 
the leading order potential, viir) for a more complicated metric, such as ([s]). Note that for 
such a metric, due to the larger number of perturbations, ^^(r) is not unique, but the equation 
Det = has several solutions v\^{r), each leading to independent linear perturbations. If any of 
f{r)v\^{r) < at some r = tq, then the static solution is unstable. We will apply this method 
to the study of instabilities for black branes in a subsequent publication. 

Now consider the case when the exact static solution, /(r), is not known. Then, finding the 
the leading order potential, VL{r), as a function of the metric components does not yet allow 
us to find instabilities. However, we can calculate VL{r) in the neighborhood of the horizon 



by making use of the horizon expansion (18). viir) can readily be obtained from (26) if we 
substitute 

Coir) -^ r^ + 24a2r2 + 72a3, 

c'o(r) -^ 4[r3-6/'(r)(a2r2-6a3)] 

c'o'(r) ^ 12{r2 + 4a2- Bras /(r) + 12a3[/'(r)]2-2/"(r)[r2a2 + 603]}, (29) 

where, for simplicity, we dropped the subscript h from r/^. We can use the expressions f'{rh) 



and f"{rh) in (29) as they are given in (18). In fact, as shown by (27), the solution is unstable 
near the horizon if the quantity f{r)vL(r) ~ (r — rh)f'(rh)vL{rh) < 0. Since r — rh > the 
solution is unstable in the immediate neighborhood of the horizon if VL{rh) = f'{Tfi)vL{rh) < 0. 
Consequently, from here on we will calculate the "effective potential" , Viirh), rather than VL{rh). 
For a2 7^ we obtain 

VL{r) = —^ ^^^ |5i?^ + 72i?^(6-43/3) 

5/R(ii2 ± 12 /? + 144 - 72/3)(i?2 ± 24 i? + 72 ^9)3 I 

+5184 i?^ (32 - 72/3 + 3/3^ - 1119744(2 - /S)/?^ ± [l728i?^(16 - 33/3) 
- 62208 /3 i?(16 - 3/3) - 60i?^] } , (30) 

where we introduced the dimensionless variables R = Th'l\a-2\ and /3 = a-ij a'l'^ and where 
±1 = Sign(a2)- We will discuss the potential at a^ = 0, which is much simpler, later. Then, as 
we pointed out in the previous section, a sufficient condition for the existence of high harmonics 
scalar instability at any given value of a^ and 03 is the existence of Re > Rmin, where Rmin 
is the minimum value of the radius of horizon, associated with the minimum mass, such that 
Vl{Rc) < 0. We will concentrate on the stability of the solution near the horizon. Note that 
inert black holes would have Rmin = 0. 
The asymptotic expansion of Vl{R), 

^^(^) ^ ]?^ > (31) 

shows that the potential is positive at large values of R. That, of course, follows from the 
stability of the Schwarzschild-Tangherlini black hole, approached by any Lovelock black hole 
when i? — )• 00. To find instabilities one needs to investigate positive, 0, and negative values of 
02 separately. 
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Figure 1. Vl{R) as a function oi R = rh' / |a2| for 02 > at two typical values of /3 = 03/02^ 



/3 = 1 (blue curve) and /3 
poleof yL(fl), Rp. 



30 (red curve). The red vertical line signifies the location of the 



As we have shown earlier, when 02 > and 03 > (20) does not forbid the radius of 
horizon, R, to approach zero. However, there is always a single positive zero, Re > 0, such 
that Vl{Rc) = and the solution becomes classically unstable at i? < Re- Re is given as 
the solution of a 6th order equation. Furthermore, Vl{R) also has a pole at < i?p < -R, 



Rp = 6(-v/2/3 — 3 — 1), if /3 > 2. The pole corresponds to the zero of CQ(r) in (26). Then it also 



has a second zero at R2 < Rp. For small R it is always negative (provided /3 > 0) 



Vl{R) 



5\/Sa3 



+ .. <0 



(32) 



Then for a^ > 2 the potential is negative at R < R2 and at Rp < R < Re- Though the potential 
is positive at the horizon when R2 < R < Rp one can show with further calculations that it 
has a range of instability above the horizon. Thus, the static solution is unstable at i? < Re- 
Whether or not Rmin < Re there cannot be any inert black hole for a > and 02 > 0. Of 
course, the exact solution shows that for this range of the Lovelock couplings Rmin = 0. The 
same thing can be discovered in a more general case, e.g. for a metric like ([8]) by connecting the 
radius of horizon with m by numerical integration. 

In Fig.l, the potential Vl{R) is plotted as a function oi R = r^^ / |a2| for 02 > and two 
typical values of /3 = 03 / 02^, /3 = 1 (blue curve), at which value the potential has a single zero 
and approaches —00 as i? — )• 0, and /3 = 30 (red curve). At the latter value of /?, besides the 
zero at R = Re there is a pole and a second zero of the potential. The location of the pole is 
shown by the vertical red line. The solution is unstable in a neighborhood of the horizon for all 
03 > 0. This has already been pointed out in P^. 

For 02 = and 03 > we obtain 



vy 



4 r [5 r^^ - 3096 r^as + 15, 552 r'^a^^ + 1,119, 744 ag'^ 
5 (r4-72a3)(r4 + 7203)3 



(33) 



where, for simplicity, as before, we dropped the subscript h from r/j. We can introduce the 
dimensionless variables p = ?^/i/(«3)^ and /x = 771/(03)^'^. Then the minimum of the mass is 



at ^min = 24^' ^ ~ 4.90 (which value could be found numerically, had we not known the analytic 
solution) attained at p = and the potential at the horizon simplifies to 



4 p'^/% p^ - 3096 p2 + 15, 552 p + 1, 119, 744] 
5{p-72){p + 72) 

The potential is positive at large p and has a zero at 



vup) = ^' '-'' ";7;^o^;:^:;^3 ' ' • (34) 



1032 96 

Pc = -^ \- — cos 

5 5 



1 1 / 45V10 

- tan"^ — 

3 I 859 



4.9769 (35) 



where it turns negative. When we further decrease p the sign of the potential changes to positive 
again at the pole p = 72 (r^ = (72)^/^ ~ 2.9129) and stays positive down to the minimum value 
of the horizon, /9 = 0. In other words, the solution is unstable at the horizon in the range 
Pp < p < Pc- Again, one can show that at p < pp there is a region of instability at non-vanishing 
radial distances (above the horizon). 

If 02 < we consider only 03 > 2a2- We use the same dimensionless quantities, /3 = 03/02^ 
and R = rj^/\a2\ as for the discussion of the 02 > case. The exact solution suggests that 
minimum value of m is 

mmin = 24a3 - 36a2, (36) 

and R = Rmin = 6 at the minimum mass. Of course, without the exact solution one can only 



find Rmin from numerical integration of ([9|). The effective potential (16) has a single pole at 

R^ -12R+ 144 - 72/3 = (37) 

The positive root of this equation is i?p = 6(1 + \/2/3 — 3) > Rmin = 6. The effective potential 
changes sign at this point. In other words, the solution becomes unstable above the horizon. 
In fact, the solution becomes unstable at an even larger value of R, because the numerator of 



Vl{R), (30) has a zero above the pole for all /3 > 2. 

Fig. 2 shows the effective potential Vl{R) as a function of R at three different choices of /3, 
the "critical" value j5 = 2 and /3 = 2.2. 2.4. In each case the potential turns negative at some 
i? > 6. In other words there is an instability above the horizon. If the black hole decays by 
Hawking radiation its radius will never reach the minimal value Rmin = 6. In the critical case 
the potential has no zero but has a pole at i? = 12. In the plot the vertical lines show the 
location of the pole at different values of /3. Summarizing our results for 03 > we found that 
in all physical cases there is a high harmonics instability at a critical mass that is larger than 
the minimal mass, whether or not the minimum radius of horizon can reach r/^ = 0. 

Finally we briefly discuss the case 03 < 0. There are three solutions of the equations of 
motion. One is asymptotically Minkowski while the other two are asymptotically AdS and dS 



(see (17) ). At a given r^ only one of them can have a horizon. Now CQ(r/j) does not have a zero 
at any Vh > 0, The potential at the horizon, Vl{R), where R = rh^ / |a2|i is positive at R > Rs 
turns negative at R = Rs where c(rs)=0 and r^ is given in ( [20| ). Vl{R) has a triple pole at that 
point. Consider first the case of Minkowski asymptotics. We proved in Sec. 3. that rh > r^. 
Vl{R) is positive for R > Rs and the solution is stable. As the radius horizon cannot approach 
zero, no inert black holes exist. Consider now case of dS or AdS. Then we must have r/j < r^. 
When r/j — ;> 

Thus at r/i — )• the solution is unstable at the horizon for 02 < 0, but stable for 02 > 0. In fact, 
in the latter case the solution is stable everywhere at small r^. However, the minimum radius 



of horizon is greater than zero. This is obvious from the exact solution, (16), and can be found 
out by integrating the equation of motion numerically. 




Figure 2. The 'effective potential,' Vl{R) as a function R = Vh' I\a2\ at the critical /3 = 2 
(red), /? = 2.2 (black), and /3 = 2.4 (blue). The location of the pole of V;,(i?), i?p, appears as a 
vertical line. At each /3 > 2 the potential has a zero at i?c > ^p- 



4. Summary 

We have introduced a simple method for the investigations of high harmonics scalar instabilities 
of static solutions in Lovelock gravity. The method is based on the truncation of the system of 
equations for linear perturbations, keeping only terms of leading order in the generalized angular 
momentum. This leads to a system of algebraic equations for the perturbations. Solving this 
system provides the leading order contribution to the potential, v^, in the Schrodinger-type 
equation satisfied by perturbations, vi^ can be calculated in the neighborhood of the horizon, 
r/i, without having an exact static solution, by applying a horizon expansion to the equations 
of motion. The terms of the horizon expansion, and consequently the leading order term of 
the potential, depend only on the radius of horizon and the Lovelock couplings. As is known 
from the investigation of Dotti and Gaiser [9| [T2] that if w^ turns negative for some r > rh then 
the static solution becomes unstable. The advantage of our method is that it can be equally 
applied to spaces with compact coordinates, in which static solutions are not generally known 
in Lovelock gravity. An example for a metric in such spaces is given in Q. 

To demonstrate the efficiency of our method we investigated the instability of static D=7 
Lovelock black holes. Since the exact black hole solution can be found jUj, the results for the 
instability of our method can be verified. In D = 7 the second and third order Lovelock terms 
contribute to the equation of motion. Both static solutions and perturbations depend on the 
two Lovelock couplings, 02 and 03. Using the horizon expansion one finds that for every choice 
of the horizon radius the solution is unique. However, from asymptotic expansions one finds 
a unique solution for 03 > 3a2^/2 only. For 03 < 3a2^ there are three asymptotic solutions, 
whereas only one of them has Minkowski asymptotics, the others are asymptotically AdS, or dS 



(see ( 17) ). However, due to the uniqueness of the horizon expansion only one of these solutions 
can have a horizon, the others must end at low r in naked singularities. These results are in 
agreement with the analysis of |17) . 

Note that our method allows one to find the instabilities as a function of the radius of horizon. 
To find the onset of instabilities as a function of the ADM mass, in cases when exact solutions 
are not known, one needs to integrate the equation of motion to find the unique relationship 
rh = F{Madm) numerically. 



We found that for 03 > the solution always becomes unstable at some Mc > Afmim 
where Mmin is the minimum mass at the given values of 03 and 02- Since the critical mass 
is determined by a high order algebraic equation, we determined the value of the critical mass 
numerically. At 03 < only a solution with a dS asymptotics has a horizon. The minimum 
radius of horizon is larger than zero. There is no high angular momentum scalar instability. 
However, the black hole with the minimal mass has a nonzero Hawking temperature and is 
consequently thermodynamically unstable. 

We conjecture that in Lovelock gravity, for all static solutions with a horizon (1) either the 
minimum radius of horizon, associated with nonzero minimum mass, (Mmin > 0) is zero, in 
which case there is a critical mass, Mc, such that all solutions Mmin < Madm < Mc have 
classical instabilities, or (2) the minimum radius of horizon is not zero and the static solution 
with Madm -^ Mmin is thermodynamically unstable. In the latter case the limiting state 
with the Madm = Mmin niust have a naked singularity. Even in case (1) a state, like the 
5-dimensional Gauss-Bonnet black hole, with r/j = would represent a naked singularity, but 
had the instability not existed, this state would have been reached by Hawking radiation only 
after an infinite time. In either case the existence of an inert massive object (with zero Hawking 
temperature) is dynamically forbidden. 

The method we use to find high angular momentum instabilities can be applied, with some 
modifications, to more complicated geometries in which one cannot even find the analytic form 
of the static solution, let alone solving the differential equation for linear perturbations, such as 
black branes in Lovelock gravity [21] [221 [23J . We intend to discuss the high angular momentum 
scalar instability of black branes ([8]) in Lovelock gravity in a future publication. Such an 
investigation may have a phenomenological interest, as well, because if in TeV scale gravity 
the radius of the horizon of a black hole is in some dimensions larger than the compactification 
radius, then the black hole appears as a black brane. 
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Appendix 

The complete contribution of linear perturbations to the equations of motion (with overall angle 
dependent factors removed) are as follows: 
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